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The Schwarzschild solution insertion in an expanding universe model, the so called "Swiss cheese 
model," is shown to possess a very unphysical property. Specifically, in this model some of the 
trajectories are discontinuous functions of their initial conditions. An alternate metric is proposed 
as a remedy. It goes smoothly between the Schwarzschild exterior solution and the Friedmann- 
Lemaitre, expanding universe metric. It is further shown that the effects of the expansion on 
planetary motions in the solar system are too small to be currently observed for this alternate 
metric. 
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I. INTRODUCTION AND SUMMARY 

Currently there are two general relativistic descriptions 
of spacetime in popular use. For planetary systems and 
other gravitationally bound structures which are small 
on the scale of the universe, there is a static description 
of the behavior of spacetime. On the other hand, for 
large-scale behavior, there is a time dependent descrip- 
tion which is appropriate as a description of phenomena 
such as the observed red-shift of distant galaxies. 

The classic question is, "How can these two disparate 
descriptions of spacetime possibly be reconciled with each 
other?" The current standard answer is that these two 
are meshed together j| on a spherical surface surround- 
ing a mass concentration which grows with time in the 
static metric, but stays at a fixed coordinate radius in 
the non-static metric. This is the "Schwarzschild solu- 
tion in a cosmological model" picture, or it is also called 
the "Swiss cheese model." This later name refers to the 
fact the in this picture, the background material is re- 
moved inside the spherical boundary. The mass removed 
depends on the mass of the central concentration and the 
curvature of space. The replacement of the mass interior 
to the sphere by a concentration of mass at the center 
is based on Birkhoff's theorem § which says that in a 
homogeneous, zero pressure cosmological model, as long 
as the material inside a sphere is spherically symmetric, 
we can replace it with a compact mass at the center with 
no change on its exterior effects. To make this matching 
work, the pressure of the exterior solution (Friedmann- 
Lemaitre metric class of the general Robertson- Walker 
line elements) must vanish, as does the pressure of the in- 
terior (the exterior Schwarzschild metric) solution. This 
condition places a restriction on the form of the univer- 
sal expansion factor of the overall universe. Although 
the matching conditions can be met, as we shall see in 
the third section, an additional problem arises when the 
dynamics are considered. 

In the second section, for the convenience of the reader, 
I gather together the necessary classical equations from 
general relativity for the study at hand. One non- 
classical item in this section is, instead of the usual 3 + 1 



spacetime split into 3 space and 1 time dimensions, I split 
spacetime into one radial coordinate, and time plus the 
two angles of spherical coordinates. The change allows 
the direct computation of the relevant extrinsic curva- 
tures. 

In the third section, I compute the stress-energy 
tensors and the extrinsic curvatures for both the 
Schwarzschild and the Friedman-Lemaitre metrics used 
in the "Swiss cheese model." By a reparameterization 
of the Schwarzschild metric, both the intrinsic and the 
extrinsic curvature can be made to be continuous. I re- 
mark that the stress-energy tensor, for certain parameter 
choices, displays no pressure discontinuity but only a cos- 
mic fluid density discontinuity. That discontinuity is in 
line with the "Swiss cheese model" idea that there are 
holes in the cosmic fluid. 

In the fourth section, I show that it can happen, for 
trajectories which approach the metric interface an near 
grazing angles, that the subsequent trajectories are dis- 
continuous functions of their initial conditions. Those 
which enter the inner or Schwarzschild metric region can 
be bound in a, finite sized, closed orbit, while those which 
do not, travel on a parabolic trajectory to infinity. To 
emphasis, this case is not the same as in Newtonian or- 
bit theory where ellipses of progressively larger size blend 
into parabolas, but here the ellipse is just finite in size! 

In the fifth section, I introduce an alternative met- 
ric. This metric is basically an adaptation of the 
Schwarzschild metric in curved space. I compute the 
stress-energy tensor. It shows an isotropic pressure, and 
no mass-density flux. It is only second order in magni- 
tude in both the Hubble constant and the inverse radius 
of curvature of the universe. The same statement is true 
of the Friedmann-Lemaitre stress energy tensor. In addi- 
tion I have computed the extrinsic curvature. Both the 
stress-energy tensor and the extrinsic curvature are con- 
tinuous, outside the Schwarzschild radius of course, as 
they come from infinitely differentiable expressions. 

In the sixth section, I compute the equations of motion 
for a freely moving test particle in the alternative met- 
ric. I then transform them to a coordinate system at rest 
at the center of the mass concentration. The flat space, 
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slowly moving particle, weak gravitational field limit of 
these equations of motion are also given. The only cor- 
rection to Newton's equations of motion in this limit is 
a term proportional to the square of Hubble's constant, 
H . 

In the final section, I gives some examples of the dy- 
namics found using my alternative metric. 



II. METRICS 

It is useful to review the properties of a few metrics. 
The line elements will all be of the general form, 

ds 2 = ^[(dx 1 ) 2 + (dx 2 ) 2 + {dx 3 ) 2 } + e^dx 4 ) 2 
= g lj dx l dx j , (2.1) 

where the Einstein summation convention is used, and 
where 

r 2 = (x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 , n = n(r,t), v = v{r,t), 

(2.2) 

appropriate to the non-static, spherically symmetric case. 
Eq. 2A defines the metric tensor gij . It will be of interest 
to know what Einstein field equation is satisfied for each 
of these metrics. The equation is, 



R 



\Rg>+kg{ = -UTi 



(2.3) 



where is the stress-energy tensor, and Rij is the Ricci 
tensor, which is a contraction of Riemann's four index 
curvature tensor. The Ricci tensor can be expressed in 
terms of the three index Christoffel symbols T as, 



Rkrn — V^kmJA (Fki),m + ^ni^km 



rm 11 1". 



hi • 



(2.4) 



where the notation )..; means take the partial derivative 
with respect to x l . In turn, the Christoffel symbols are 
defined in terms of the metric tensor as, 

T™ = \g mk { 9 M,3 + 9kj.i - 9ij,k), where g mk g kj = SJ\ 

(2.5) 

defines g km , and 6™ is the Kronecker delta function. Fi- 
nally, 



(2.6) 



is the contraction of the Ricci tensor. 

The T 44 element is the mass-energy density. The T 4/3 
is the mass-flux density through an area perpendicular 
to the direction (3 per unit time. (Greek indices run from 
1 to 3 while Roman indices run from 1 to 4.) The T aa 
element is the pressure in the a direction, and the T a/3 
element is the flux density of the a component of mo- 
mentum in the (3 direction. It is manifest that eq. 2.3 
allows the direct computation of the stress energy tensor 



from the metric tensor. Tolman || gives the results for 
the form of the line element 



ds 2 



-e^(dr 2 +r 2 d8 2 + r 2 sin 2 



e v dt 2 , 



(2.7) 



which is the spherical coordinate version or eq. 2.1. The 



non- vanishing elements of the stress-energy tensor are 
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fj, 



where A is the cosmological constant, an over dot denotes 
the time derivative, and a prime denotes the derivative 
with respect to r. 

Another quantity which is useful to consider is the ex- 
trinsic curvature. Q This concept arises in the Arnowitt 
et al. II splitting of four dimensional spacetime into 3 
dimensional space plus one dimensional time. The rela- 
tionship between the metrics is given by 

ds 2 = ■ g,, ,!.r' d.r' (2.9) 
= {3) g a f3{dx a + N a dt)(dx + N f3 dt) + N 2 dt 2 , 

where the N a are the three shift functions and N is the 
lapse (of proper time) fu nctio n. The intrinsic curvature 
is the analogue of R [eq. 2.6] in three dimensions. The 
extrinsic curvature measures the fractional shrinkage and 
deformation as one advances in time from one space-like 
hyperplane to the next. The extrinsic curvature tensor 
is given as 



I 

2N 



N, 



a\/3 



dg a 



Of 



(2.10) 



where the notation )i a means the covariant derivative 
with respect to x a . 

The reason for the interest in the extrinsic curvature in 
our case is that the necessary and sufficient junction con- 
ditions || to join two metrics in 4 dimensional spacetime 
across a 3 dimensional hypersurface is that both g a p and 
K a p be continuous across the surface. In our ca se w e are 
concerned with line elements of the class of eq. [2.7| . The 
split is between r instead of t and the other three vari- 
ables. In this case the shift functions are all zero, and the 
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lag function is N 
K are 



exp(0.5/i). The non-zero elements of 



K- 



22 = i(A 2 + 2r)exp(0», 
K 33 = i(/i'r 2 + 2r)sin 2 6'exp(0.5^), 
K44 = —\v' exp(^ — 0.5/x). 



(2.11) 



Returning to the three-space, constant-time formalism, 
we give the equations of motion of a free test particle as 
seen by local co-moving observers along the test particle's 
path. That is to say, a set of observers whose coordinates 
do not change with time. In other words, we need the 
equations for the geodesic curves in spacetime. For the 
class of line elements we are considering, it is sim plest to 
start with a Lagrangian formulation. By eq. 2.1 we may 
write this formulation as, 



= / Ldt,=> L = — = (gij: 



,1/2 



(2.12) 



The standard Euler-Lagrange equations for an extreme in 
path length (Here we seek a minimum distance between 
two fixed endpoints.) are 



d ( dL\ _ dL 
di \ ~d±" / ~ ~dx" 



(2.13) 



Thus, using the diagonal nature of the metric tensor, we 
obtain for the standard geodesic equations, 



ds\ d 
dt ) di 



where, of course dx /dt 
get 



1 dx 1 dx^ 

2-'"-^r^r- (2 - 14) 



1. From the line element we 



|) = [e"+g a0 i a ±P] 1/2 



(2.15) 



These results display, for the class of line elements we are 
considering, the three, second-order, non-linear, coupled 
equations for the three coordinates x a of a test particle 
as a function of time. We have assumed isosynchronous 
coordinates everywhere in the three-dimensional, space- 
like hyp er-surface. The square roots can removed from 
eq. 2.14 by rewriting it as 



ds 
dt 

d 
di 



d_ 

dt 



[g ai x l ] 



1 dx 1 dx^ 



(2.16) 



III. "SWISS CHEESE MODEL" 

As mentioned in the first section, there is a popular 
cosmological model which in the large has the Friedmann- 
Lcmaitre line element, 



ds 2 = - 



a{tf 



„3\2i 



[1 + (r/2i?) 2 ] 2 
c 2 (dx 4 ) 2 



[(dx'Y + {dx 2 ) 2 + (dx J ) 



(3.1) 



where in the notation of Sec. II, 

«W 2 



[1 + {a{t)r/2a{t)R) 



2' 2 



(3.2) 



where c is the velocity of light and a{t)R is the radius of 
curvature of the model universe. 

The non-zero elements of the stress-ene rgy tensor as- 
sociated with this line element are, by eq. 2.8, 



8irTl = 8ttT 2 2 = 8ttT 3 3 = 



1 



[a(t)R\< 



[a(t)R}< 



-A 



(3.3) 
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A = 87T/900 



The extrinsic curvature tensor, as given by eq. 2.11 for 
this metric has the non-zero elements, 



#22 = 



a{t)r [l - (r/2Rf 



[1 + (r/2R) 



2] 2 



K 33 = sin 2 6K 



22- 



(3.4) 



In this cosmological model on scales small compared 
to that of the universe, as in, for example, the solar sys- 
tem, a quite different metric is used. It is Schwarzschild's 
exterior solution. The line element for it is 



ds z 



(l+£) [(^) 2 + (^ 2 ) 2 + (dx 3 ) 2 ] 



1 — m/2r 



1 + m/2r 
where in the notation of Sec. II 



{dx A )\ 



(3.5) 



^ = 1 1 + — 
2r 



1 — m/2r 
1 + m/2r 



(3.6) 



Here m is an abbreviation for GM/c 2 where M is the 
mass concentration. This metric is only valid outside the 
Schwarzschild radius r >> r$ — 2GM/c 2 . Inside this 
radius, Schwarzschild's interior solution is required but 
we shall not be concerned with this aspect here. 

There are no non-zero elements of the stress-energy 
tensor associated with this line element, unless the cos- 
mological constant is different from zero. Thus it cor- 
responds to zero pressure and zero density, except for a 
mass concentration in the center. If A =/= then, 



T l = T 2 = T 3 = T t = -A/87T 



(3.7) 
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This result would correspond to a uniform density and 
pressure through out space, rather than the empty space 
with A = 0. 



A cross comparison of eq. 3.7 with eq. 3.3 shows that 
unless, 



1 



[a(t)R}< 



ac I ac z 



(3.8) 



the stress-energy tensor is discontinuous at the boundary. 
The only solution of these equations, as R is a constant, 
is a is a constant, i.e., a linear expansion factor, or a 
static flat universe. 

Before we can us eq. 2.11 to compute extrinsic cur- 
vature for the Schwarzschild case, we must first repa- 
rameterize the metric so that the "Swiss cheese model" 
metric boundary is given by one of the coordinates equals 
a constant. Q It will be convenient to use spherical co- 
ordinates. If we follow the textbook [Q approach, then 
we should start with a zero pressure cosmological model. 
By Birkoff 's theorem in this case we may hollow out a 
sphere and replace it by a mass concentration at the cen- 
ter. Since the density of the "cosmic fluid" is inversely 
proportional to a 3 , the radius of the sphere is fixed in 
the Friedmann-Lemartre coordinates. The surface of the 
sphere is a three dimensional hypersurface parameterized 
by the time and the two angle variables of spherical coor- 
dinates. To proceed, we note that an observer sitting on 
the boundary is on a geodesic for rp^ = a constant for 
all time. The relation between the Friedmann Lemaitre 
time i and the Schwarzschild variables as seen by this 
observer is 



{dsf = c 2 {dt e f 



1 - 



2r K 



1 



m 
2rT 



- 1 



m 
2t> 



dt P 



(3.9) 

{dt e )\ 



Since the observer is on a geodesic, we m ay deduce, by 
means of the time component of eq. 2.14, that he sees 



1- — \ 

2r 1 & _ K i/2 



(3.10) 



where if is a constant of integration. Note is taken that, 
as the integration is ov er t, K may d epend on r, of course. 
By combining eq. 3.9 and eq. 3.1C, we obtain, 
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which gives the behavior of r e , t e as a function of the 
Friedmann-Lemitre time variable, as seen by the comov- 
ing observer sitting on the metric boundary in the Swiss- 
cheese model. 

By equating the coefficients of the angular variables, 
we obtain, 



m 

2r~ e 



a(t e )r e 

+ (— Y 

\2RJ 



(3.12) 



d(t e ) 
a(i e ) 



1 dr e 

r R dt e . 



1 - 



2r fi 



777 

2rT 



dtg 

dL 



Note is taken that the boundary coordinate f e in the 
Fricdmann-Lemaitre metric depends on the central mass 
concentration, and on the radius of curvature of the uni- 
verse, i.e. on 777 and R. 



Next we combine eqs. 3.11 and 3.12 and comparing 



the result with the "cosmological equation." This equa- 
tion is the equation for the X4 component of the stress 
energy tensor derived from the Friedmann-Lemaitre line 
element. The result of the application of this textbook 
method is that the constant of integration is 



K 



(3.13) 



which is a function of m and R. 

We are now in a position to introduce a reparameter- 
ization of the Schwarzschild metric. We choose the new 
parameters t = t and f. The second variable is chosen 
so as to keep the metric continuous at the interface, and 
to keep the metric diagonal, if possible. Th e re para mc- 
terized Schwarzschild metric is given by eq. 3.1 , Eq. 3J 



insures that 344 = c . The equation of continuity, and the 
vanishing of the elements gi4 = gn yield the conditions, 



dr, 
df, 
dr\ dr 
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(3.14) 



The solution of these equations is 



dr 
df 



K l / 2 a{t) 
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Arl 
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dt 
df 



K^ 2 a(t) ( 1 



m \ 3 
2r, 



V2i?/ 



l- — 



i - 



4r? 



1 - 




1/2 



(3.15) 



3ii 



2 \ 2 
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1 - 



4r 2 



V2i?/ 
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Our construction insures the continuity of 511 = g\\ at 
the metric interface. The conti nuity of the (22) and the 
(33) elements are insured by eq. 3.12. All the off-diagonal 
elements vanish. 



We may now apply eq. 2.11 to o btain the extrinsic cur- 
vature. The result is just eq. 3.4 with (r, t, 9) replaced 
by (f, t, 9). This result shows the continuity of the 
extrinsic curvature. 

There is however, one important item to note. Con- 
sider the derivatives perpendicular to the boundary hy- 
persurface, evaluated at the interface. 

dlog(-ffu) r 



dr 

dlog(-gn) 
df 



i? 2 + 0.25f 2 ' 



«(*): 



i-f-V 

\2RJ 



m 
4r 2 



(3.16) 



To leading order in a flat universe, eq. 3.16 becomes, 



dlog(-gn) = 
df 



rflog(-gn) 
df 



(3.17) 



That such a discontinuity may occur is well known [(7) . I 
will explore some of the consequences of this discontinu- 
ity in the next section. 

As was remarked above, it is well known that it is a 
necessary condition that the pressure equal zero in order 
for this matching to occur pj. As the pressure p = —T ( c * 
for each a (not sum med here), the zero-pressure con- 
dition means, by eq. 3.7 which gives the pressure for 
the Schwarzschild case in terms of the T£ elements that 
A = 0. Turn ing to the the Friedmann-Lemaitre case, we 
see from eq. 3.3 that it must be that 



1 



[a(t)R}< 



= 0. 



(3.18) 



(t/t 



If use su bstitute the standard form a(t) 
eq. 3.18, then the only solutions are R — 00 together 
with ip — or 2/3. The first solution is the trivial case 
of a static universe and is of no interest in the present 



discussion. In the second case, the stress-energy tensor 
element T4 > 0, which means there is a mass discontinu- 
ity. This is expected by the structure of the approxima- 
tion of scooping out a hollow sphere, and then having a 
uniform density outside. 



IV. DYNAMICS IN THE "SWISS CHEESE 
MODEL" 

Following up on the discontinuity in the derivate of 
the metric, perpendicular to the metric boundary in the 
"Swiss Cheese Model," we investigate the some of the 
dynamic properties of this model. It suffices for my pur- 
poses to consider the simpler case of a flat (R — 00) uni- 
verse. The equations of motion are as follow s. For the 
Schwarzschild metric, we get from eq. 2.14 , to leading 
order, the well known Newtonian equations 



cPp 
dt 2 



GM _ 



(4.1) 



The solutions to this equation are the familiar Newtonian 
conic sections. 



For the Friedmann-Lemaitre metric from eq. 2.14 we 
get a first intergal as 



A 



(4.2) 



For my purposes, it is more convenient to use a vari- 
able more closely equal to the proper distances. Thus for 
p = a(t)x, eq. 4.2 becomes, 



dt 2 



a(t) 
a(tY 



(4.3) 



In both metrics p is within plotting accuracy for the 
proper distance in the examp les I will consider. 
The general solution of eq. \L3 is 



At* + Bt 1 ^ = At 2/3 + Bt 1 / 3 , 



(4.4) 



when the currently, theoretically favored value of ip = | 
is chosen. Let us take the example in rectangular coordi- 
nates where initi ally x = A, x = 0, y = 0, y = X at time 
t = to. Then eq. 4.4 becomes, 



A 



y = 3£ A 



2/3 



2/3 



2 - 



1/3 



1/3' 



(4.5) 



It is worth noting that for t — to <C to that the motion is 
almost exactly that of a straight line, as is to be expected 
in flat empty space. Specifically, direct calculation yields 



S(t ) 



2A 

'9tf' 



y(t ) = 0. 



(4.6) 
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This apparent acceleration is of order Hq . 



The solutions in eq. 4.4 are quadratic in the parameter 
t 1 / 3 . Thus we can obtain the equation for the trajectory 
as a quadratic plus linear expression in p x and p y . The 
form would be 



c(ap x + bp y ) 2 + d(ap x + bp y ) 



(4.7) 



This form is readily recognized as a parabola. Thus it is 
the case that a freely moving particle in an Friedmann- 
Lemaitre expanding space always appears to be moving 
in a parabola. This effect is c ause d by the small d term 
which appears in the equation 4.3 as a forcing term. 



As an illustration of the behavior of the "Swiss cheese 
model" I have computed the following trajectories. I use 
as a unit of time the Hubble time, that is 1/Hq, which 
is of t he order o f 10 18 seconds. As a unit of distance I 
use \J Mq G I Hq which is about 20 million Astronomical 
units. Mq is the mass of the sun. I set to — 1 to switch 
to our current units. The metric interface is at r = t 2 ^ 3 
in these units for flat spacetime, as mentioned above. In 
order to follow a Friedmann-Lemaitre trajectory the test 
particle's distance from the origin must be larger at every 
time than that for the interface. Thus, 



10t 4/3 - 22i + 13t 2/3 > t i/3 



(4.8) 



The trajectory will intersect the interface if eq. 4.8 is an 
equality. By means of the quadratic formula, an inter- 
section will occur if 



t l/3 



11 ± v/-9 + 13/A 2 
10- 1/A 2 



(4.9) 



It will be observed that for A < V13/3 there are two 
real roots. If A = 1, then t 1 ' 3 = 1, 13/9. On the other 
hand, if A > VT3/3, the roots are imaginary, so there are 
no intersections. If A = VT3/3 there is a double root at 
t 1 ' 3 = 13/11. In this case the parabolic trajectory just 
grazes the metric interface. 



Friedmann-Lemaitre ' 




point of separation v 



-5 

X 



FIG. 1 "Swiss cheese model" trajectories which begin 
at xo = \/l3/3 + e (the Friedman-Lemaitre case), and 
begin at xo = VT3/3 — e (the Schwarzschild case). Here, 
e > may be chosen as small as one pleases. The ini- 
tial part of both trajectories is a parabola generated by 
the Friedmann-Lemaitre metric. At time t = (13/11) 3 , 
marked in the figure, the two trajectories separate. 

In Fig. 1, 1 illustrate the two different trajectories when 
ro = VT3/3 and 0o = 1-0. In the case where ro is just any 
arbitrary amount smaller, the expanding interface over- 
takes the test particle following its Friedmann-Lemaitre 
parabolic trajectory and it must then follow the static 
Schwarzschild equations of motion. The Schwarzschild 
metric takes over at t = (13/11) 3 as explained above, 
and after that the trajectory is an ellipse with scmimajor 
axis 7.6287 and the semiminor axis 3.986. These imply 
an eccentricity of 0.69068 and a semilatus rectum pa- 
rameter of p = 2.086. On the other hand, if ro is any 
arbitrary amount larger, it escapes the moving interface 
and continues to follow the parabolic trajectory. It is 
evident from Fig. 1 that future trajectories are, in some 
cases, discontinuous functions of the initial conditions for 
the "Swiss cheese model." 

Put another way, the Schwarzschild metric permits 
closed orbits and the Friedmann-Lemaitre metric does 
not. In terms of the latter coordinates, one can choose 
initial conditions so that the parabolic trajectory just 
grazes the metric interface (fixed radial coordinate in this 
metric) and the speed is low enough that for infinitesi- 
mally different initial conditions the trajectory crosses 
the interface and is caught in a bound state, or alterna- 
tively misses the interface and proceeds on its parabolic 
trajectory. All these effects take place in supposedly 
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empty space of the order of 20 million AU from a mass 
concentration of size M Q , and are quite counter to one's 
physical intuition that such discontinuities should not oc- 
cur there. 



V. AN ACCEPTABLE METRIC 

I propose the following line element to represent a mass 
condensatio n in an expanding and curved universe. It is 
of the form 2.7 where 



a(tf 



[1 + (a(t)r/2a(t)R) 

1 - m/2a(t)r 
l+m/2a{t)r 



212 



2a(t)r 
GM 



1 



(—)' 

\2RJ 



(5.1) 

n 1/2 



where G is Newton's constant of gravitation. It is 
to be noted that this metric is an adaptation of the 
Schwarzschild metric in curved space. It is not claimed 
that this metric is unique. Certainly any coordinate 
transformation of this metric is equivalent. It does how- 
ever have the property that in the limit where a(t) is 
a constant, it reduces to the Schwarzschild metric in 
curved space. Also, when the central mass vanishes it 
reduces to the Friedmann-Lemaitre metric. It is an in- 
finitely diffcrcntiable solutions to the Einstein field equa- 
tions which corresponded to values of the stress-energy 
tensor which are only of second order in and Hq, 
except at the central mass. This latter property is also 
true of the Friedmann-Lemaitre metric in curved space, 



as may be seen in eq. 3.3. Thus we have, by this example 



whose properties are globally similar to the Friedmann- 
Lemaitre metric, demonstrated that the Swiss cheese 
model is not require to match the expansion of the uni- 
verse observed at large scales and the absence of any such 
expansion observable in the solar system. As shown in 
the previous section since the Swiss cheese model has a 
rather severe and very unphysical deficiency, we think it 
is time to begin the search for an acceptable metric. 

The non-zero elements of the extrinsic curvature is 
given by (2.11) as 



K 33 = sin 2 9K 22 

mc 2 (l — m/2a(t)r) 



(5.2) 



K 44 = — 



a(t) 2 r 2 (l + m/2a(t)r) 5 

In the lim it that m — > these curvatures reduce to those 
of eq. 3.4 and in the limit R — > oo and a = (a = 1) 



they reduce to 



K 2 2 = r 



l-(™ 

2r 



K 33 = sin 2 6K; 



22- 



K i4 = - 



mc 2 (l — m/2r) 
r 2 (l + m/2r) 5 



(5.3) 



These curvatures agree with those of the unreparame- 



terized Schwarzschild metric. The curvatures in eq. 5.2 
differ only in that r is replaced by a(t)r, R by a(t)R, and 
there are corrections for the overall curvature of space. 
Thus, it reflects, as does the metric, the very same be- 
havior, in terms of a(t)r instead of r as was found by 
Schwarzschild for his metric. 

Both the metric and the extrinsic curvature are contin- 
uous outside the Schwarzschild radius, which is necessary 
for a metric to be acceptable. 

We now turn to the stress-energy tensor. The current 
me tric shares with both the Friedmann-Le mai tre metric, 
eq. 3.1, and the Schwarzschild metric, eq. 3.5, the prop- 
erty that Tf = T\ = 0, so there is no mass-density flux 



This pr oper ty follows by direct computation from eq. 2.8 



and eq. 5.1. For the other non-zero components, we find 



1 + {m/2a(t)r)f 



a{t) 2 R 2 {l + m/2ar) 5 (l - m/2a{t)r) 
2fia[l + m/2a{t)r)] + 3d 2 [l - m/2a(t)r] 



+ 



a{t) 2 c 2 {\ - m/2a(t)r) 
8ttT 2 2 = 8ttT 3 3 = 8nT^ 

3 3d 2 



-A 
(5.4) 



8nT? 



a{t) 2 R 2 {l + m/2a{t)rf 



A, 



In the limit as m — > o r r — > oo these tensor elements 
reduce to those of eq. and also in the limit as R — ► oo 
and a — > they vanish as with the Schwarzschild metric, 
unless A 7^ 0, in which case the result (3.7) is obtained. 
The dominate terms in the diagonal elements (pressure 
and density) are a sum of terms of second order in 1/ R 
and terms containing two time derivatives of the univer- 
sal expansion factor a(t). The metric given by eqs. 
and 2.7 is the solution of Einstein's field equations eq. 



when the stress-energy tensor, eq. 5.4, is specified. 



5_1 
2.3 



VI. DYNAMICS 

In this section we investigate the equations of motion 
a t est p article i n o ur proposed metric as described by 
eq. [2.7| and eq. 5.1. Our treatment diverges from that 
of McVittie H at this point as he directly substitutes 
p = a(t)r into the line element rather than using the 
transformation equations 



9ki 



dx l dx^ 
dx k dx l ' 



(6.1) 



The correct change of the line element for this change of 
variables is 

ds 2 = -e u [(dp) 2 + P 2 {d0) 2 + p 2 sin 2 6(d(j>) 2 ] 

, - -e u (^Y (dT) 2 + 2e u (^)dpdr, (6.2) 
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where 



v 2 

e = c 



(1 + m/2p) 4 
[1 + (p/2a{r)R) 2 } 2 

1 — m/2p 



1 + m/2p 



(6.3) 



Clearly these are non-synchronous coordinates as the co- 
efficient of dpdr is non-zero. 

We will make this change of variables later after the 
equations of motion have been derived. The main differ- 
ence is that McVittie omits the p terms from his subse- 
quent equations. The equations of motion then become, 
using eq. 2.14. for 9 and </>, 



ds\ d 
'dt ) ~dt 



ds\ d 
di) dt 



dt 



e M r 2 sin 2 I 



= e'V 2 sin 9 cos t 



= 



(6.4) 



One can see by inspection, the motion in the plane 
9 = 7r/2 is a solution. The intergal of the second equation 
gives the result, 



r 2 <j) 



Ae~ 



(6.5) 



where A is a constant of integration. This equation is 
the conservation of angular momentum in this coordi- 
nate system. To obtain the equation of motion for r, it 
is convenient to use eq. 



d 
dt 



2.16j . We obtain, 

dt J dt 1 1 

p'e^r 2 + (p'r 2 + 2r)e^9 2 
+ 2r)e^ sin 2 9 ft - v't 



where, from the line element, 



r 2 sin 2 i 



(6.6) 



(6.7) 



At this point, we simplify to motion in the 9 — tt/2 



plane. Thus 9 = 9 = 0, and eq. 6.6 becomes 



ds\ d 
dt) ~dt 

1 

2 



dt 



p'r 2 + (p'r 2 + 2r)4> 2 



(6.8) 



These formulas yield the equations of motion of a test 
particle in the reference frame which is at rest with re- 
spect to the coordinate system at the location of the test 
particle. We want to find the motion with respect to an 



observer at rest at r = 0. To this end we introduce the 
change of variables, 



P 

a (ty 



p 



a{t) a(t) 2 



(6.9) 



This change of variables yields coordinates which are 
equal to the proper distances, as viewed from the ori- 
gin, when t he m ass concentratio n is absent. W ith the 
substitution 6.£ we obtain for eq. |6.5| and eq. 6.7 



/rV ' = A " (/,2 <~" ( J)' 



ds\' 
dt) 



Finally, eq. 



„<2 . 

p-2-p 
a 



a{tf 



P- 



a(t) 



+ p- 



212 



becomes, 

2[^V-- 
a I a 



Ml P~ ~P 



1 , (P 
—pa 



TP 



1 t ,P 



a a"' J 2 v a 

1 fds\ 2 d 

2 {'dt ) di 



A 2 a i e- 2 » 



(6.10) 
(6.11) 

(6.12) 



a 

P--P 
a 



where we have used eq. 6.10 to eliminate the <f> depen- 
dence. There is also a <p in ds/dt but it too can be elim- 
inated by the substitution of cj) from eq. 6.1C in eq. 3.11 . 
To assess the importance of the various terms it is helpful 
to note the following dimensionless quantities, in "plan- 
etary units," 



T@H 

GM Q 
c 2 R ffi 



5X10- 11 , ^ 



1 x 10" 



c 

( R 



bbk 



1 X 10" 

2 



0.6 x 10 



-30 



(6.13) 



where 7®, i>®, i?® are the orbital period, velocity, and 
radius of the earth, and M® is the mass of the sun. Some 
limiting cases are of interest. First, we take the flat-space, 
slow-speed limit, i.e., R — oo, and 



/ 1 GM\ 

2c 2 p 
GM 
V 1+ 2cV/ 



(6.14) 



Thus eq. 6.12 reduces to 



p- -p- yp- - f > 

m{p — ap/a 



a \ dm(3 — m/p) 



ap[l - {m/2p) 2 } 
A 2 c 2 {\-m/2pf 



p 2 (l + m/2p) p 3 (l + m/2py 



(6.15) 

mc 2 (l — m/2p) 
p 2 (l + m/2p) 7 



We may further reduce these equations by discarding 
terms which are proportional to c~ 2 for when t he v eloci- 
ties are much less than the speed of light. By eq. |5.l| these 
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are the terms proportional to m al one, but we retain, of 
course, the terms in mc 2 . Eq. 5.15 reduces further to 



2-2 



d A 2 c 
~P = 



GM GM 



P* 



P^ 



- 1 



(6.16) 



The coefficient of the discarded term on the left-hand side 



of eq. 3.15 is of the order of 10 19 as it is the product 



of two first order corrections. The discarded (first) term 
on the right-hand side of eq. 6.15 is smaller by a factor 
of p 2 /c 2 than the last term. The other items discard are 
factors of GM/c 2 p smaller than the dominant term s. 

that 



It is to be noticed that in the limit of eq. 6.16 



it differs from Newton's equation of gravitation only by 
a term proportional to d. Note is taken the Noerdlingcr 
and Petrosian ]9| do take account of this term. The con- 
stant po is just another form of the constant of integration 
A. The magnitude of the d term equals that for the sun's 
gravity at about 0.5 kiloparsecs. The effects on the scale 
of the solar systems are too small to be measured. 



VII. EXAMPLES 

The equations of motion in a flat, Friedmann-Lemaitre 
expanding universe for a slowly moving test particle un- 
der no external force s ar e, by eq. 6.16 and the correspond- 
ing reduction of eq. p.% 



dt 



p 2 =0, (7.1) 



for the standard form of the universal expansion factor 
for the universe, as describe at the end of section III. t 
is the current age of the universe. One easily recognizes 
these equations to be just exactly Newton's equations in 
a plane in spherical coordinates. If w e ch ange to rectan- 
gular coordinates, we get exactly eq. 4.3. The behavior 



of the solutions is discussed in detail in Section IV above. 

Next we consider the case where we add gravitational 
effects to the ir lea ding order. For purely radial motion, 
the A of eq. 6.10 is zero. Thus the equation of motion 



6.16 becomes, 



P 



a(t) 
o(t) 



P - 



GM 



(7.2) 



This equation differs from the Schwarzschild metric equa- 
tion, 4.1 by the addition of a term in d. We will consider 



this behavior over time periods short compared to the 
age of the universe. If we multiply by dp and integrate, 
we get, 



1 



'la 



P 



GM 
P 



E 



(7.3) 



where Eq is the constant of integration. If d = 0, then 
Eq = would correspond to a test particle which had 
zero velocity at p — oo. I choose to examine this special 
case. Then eq. 7.3 becomes 



PO 2 y/GM + %? 

" ■ -i ( e /2 x ' 

-sin = 

a \ y/-2GMa/d 




t-tn 



t-tn 



(7.4) 



/HI 



by Pierce's tables pfl. Thus, 



P 3/2 = pT + 



Pq' 2 + -V2GM 



2GM a , 3, . a 



(t-to) + p{t-t a ) 3 + 



(7.5) 



It is evident from this solution that the leading order 
corrections due to the expansion of the universe are of 
the order H^t — to) 2 which is extremely small on the 
planetary time scale. The solution when (t — to) = O(to) 
would take account of the time dependence of d/a. The 
more general case, where E 0, can also be integrated 
in terms of elliptic functions of the first and third kinds 

Next we invest igate bound circular motion. To do so, 
I set p — in eq. 5.16] which gives, 



2„2 



A z c 



GM 



(7.6) 



If I use eq. |6~10| to reintroduce <j) 1 and remember that the 
period T = 2ir /<j), then I find, 



2^3 



4ir z p 
GMT 2 



ap 

GMa 



(7.7) 



which is Kepler's law relating the square of the period to 
the cube of the radius, with a correction caused by the 
expansion of space. For the case a(t) cx t 2 ^ 3 the correc- 
tion term becomes, —^H'qP^/GM. If we use the solar 
mass, then 



2„3 



47rV 
GMT 2 



= 1 - 3.307 x 10" 23 ^ p 3 , 



where p is in astronomical units and /150 
Hq = 50 km per second per megaparsec. 



(7.8) 



1 when 
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1 ' \ 

Friedmann-Lemaitre 



Schwarzchild 




GM 



A 



i GM 



A 



FIG. 2 Trajectories in the static (Schwarzschild) met- 
ric and the non-static (Friedman-Lemaitre) portions of 
the "Swiss cheese model." The trajectory for the metric 
of this paper, with the same initial conditions as for the 
Schwarzschild case, is also shown. 

These results are illustrated in Fig. 2. In flat space at 
time to, the interface between the static Schwarzschild 
metric and the non-static Friedmann-Lemaitre metric is 
at distance unity from the origin, in the units of sec- 
tion IV. We display the Schwarzschild result when started 
with unit angular velocity (0=1). It is a circle of radius 
unity. Next we display the trajectory just outside the 
interface. It is, as expected, a parabolic curve. I have 
started it with unit angular velocity (again <j> — 1) at a 
distance of 1.25 from the mass concentration. It is to be 
noted that after a unit (Hubble) time has passed, the two 
trajectories are significantly separated. 

The reason for starting it further out is, as explained 
in section IV, that if I were to have started it at a dis- 
tance between 1.0 and V43/3 ps 1.2018504, the trajectory 
would have been over taken by the expanding spherical 
interface between the two metrics of the "Swiss cheese 
model." 

In addition I show the trajectory using the presently 
considered metric. I have again started with unit dis- 
tance and unit angular velocity. The result here is that 
it converges fairly quickly to an elliptical trajectory. In 
this case I find, using the standard equations, 



(7.9) 



where A, B are the major and minor semi-axes, and 
the areal velocity A is a constant by the conservation 
of angular momentum. In this case I find A ~ 1.0437 
and B ~ 1.0216. In as much as the semilatus rectum 
parameter p = 1, the latus rectum itself is clearly de- 
fined by the intersection of this ellipse with the unit cir- 
cle (Schwarzschild trajectory). The latus rectum is the 
line through the focus (origin in this case) which is per- 
pendicular to the semimajor axis. The eccentricity is 
e fa 0.20462. 
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